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ABSTRACT 



Context. Over the last few years, the chemistry of molecules other than CO in the planet-forming zones of disks is 
starting to be explored with Spitzer and high-resolution ground-based data. However, these studies have focused only 
on a few simple molecules. 

Aims. The aim of this study is to put observational constraints on the presence of more complex organic and sulfur- 
bearing molecules predicted to be abundant in chemical models of disks and to simulate high resolution spectra in view 
of future missions. 

Methods. High S/N Spitzer spectra at 10-30 ^m of the near edge-on disks IRS 46 and GV Tau are used to search for 
mid-infrared absorption bands of various molecules. These disks are good laboratories because absorption studies do 
not suffer from low line/continuum ratios that plague emission data. Simple LTE slab models are used to infer column 
densities and excitation temperatures for detected lines and upper limits for non-detections. 

Results. Mid-infrared bands of HCN, C2H2 and CO2 are clearly detected toward both sources. As found previously 
for IRS 46 by Lahuis et al. (2006), the HCN and C2H2 absorption arises in warm gas with excitation temperatures of 
400-700 K, whereas the CO2 absorption originates in cooler gas of ^250 K. Absolute column densities and their ratios 
are comparable for the two sources. No other absorption features are detected at the 3a level. Column density limits of 
the majority of molecules predicted to be abundant in the inner disk — C2H4, C2H6, CgHe, C3H4, C4H2, CH3, HNC, 
HC3N, CH3CN, NH3 and SO2 — are determined and compared with disk models. Simulations are also performed for 
future higher spectral resolution instruments. 

Conclusions. The inferred abundance ratios and limits with respect to C2H2 and HCN are roughly consistent with 
models of the chemistry in high temperature gas. Models of UV irradiated disk surfaces generally agree better with the 
data than pure X-ray models. The limit on NH3 /HCN implies that evaporation of NHs-containing ices is only a minor 
contributor. The inferred abundances and their limits also compare well with those found in comets, suggesting that 
part of the cometary material may derive from warm inner disk gas. The high resolution simulations show that future 
instruments on JWST, ELTs, SOFIA and SPICA can probe up to an order of magnitude lower abundance ratios and 
put important new constraints on the models, especially if pushed to high S/N ratios. 

Key words. Protoplanetary disks - Astrochemistry - ISM: molecules - Line: profiles - Planets and satellites: formation 



1. Introduction 

The chemical composition of the gas in the inner re- 
gions of circumstellar disks plays an important role in 
determining the eventual composition of the comets and 
atmospheres of an y planets that may form from that 
gas (see reviews by Primil 1993t jEhrenfreund fc CharnlevI 

23 [ 



120001: iMarkwick fc Charnlevll2004tlBerginll2011[ ). In the last 
few years, observations with the Spitzer Space Telescope 
have revealed a rich chemistry in the inner few AU of 
disks around low-mass stars, containi ng high abundances 



of HCN, C2H2, CO2, H2O and OH (iLahuis et al 



Carr fc Naiita"2008':'Salvk et al.ll2008': '"Pascucci ct al 



2ooe 



Carr fc N aiita 2011: Kr ueer et al.ll20 11: Naiita ct al. 2010; 



200c 



Pontoppidan'et al.il2010HSalvk et af.lf201llD . Soectrallv and 



spatially resolved data of CO using ground-based infrared 
telescopes at 4.7 /xm show that the warm molecular ^as is 
indeed a ssociated with the disk (e.g., Naiita et al. 2003; 
Brittain ct al. 2003. 20071 l2009t iBlj : |kc fc googcrt 20[)4: 



Pontoppidan et al.ntoodSalvk et at[l2'oilU iBrown et all 



submj), with in some cases an additional contribution from 



a disk wind (|Bast et all I2OIII: iPontoppidan et all I2OIID . 
Spectrally resolved ground-based obser vations have also 
been obtained of OH and H2O a t 3 fim fCarr et al."200j; 
MandcU ct al. 2008; Sal vk et al.l [2b08: Fc dcle ct al. 2pij|, 
and most recently of HCN and C2H9 (ICibb et al] 120071: 
iDoppmann et al.ll2008t iMandeh et alll2012[ ). Ah of these 
data testify to the presence of an active high-temperature 
chemistry in the upper layers of disks that drives the for- 
mation of OH, H2O and small organic molecules. However, 
it is currently not known whether this gas contains more 
complex organic molecules which may eventually become 
part of exoplanetary atmospheres. 

Observations of large interstellar molecules are usually 
carried out using (sub-)millimeter telescopes. A wide vari- 
ety of complex organic species have been found in low- and 
high-mass protostars at the stage when the source is still 
embedded in a dense envelope (see iHerbst fc van Dishoecld 
2009, for review). For disks, the pure rotational lines of CO, 
H2O, HCO+, H2CO, HCN, N2H+, CN, C2H, SO, DCO+ 
and DCN have been reported but more complex molecules 
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have not yet been detected (e.g., iDutrev et al.l 

Knstnnr pt nl I ll QQf- 'Th\ pt nl ' 9004'- Vimntp nt al I 



1997 



201C; 



Kastner et all llQQT": 'Thi et al.' 2004; Fuentc et al 

Hcnn i ng et al.l [2010: Oberg et al. 2011; Hoaerheiid e et al.l 
201 ih . Although these millimeter data have the advantage 



that they do not suffer from dust extinction and can thus 
probe down to the midplane, current facilities are only sen- 
sitive to the cooler gas in the outer disk (> 50 AU). Even the 
Atacama Large Millimeter/submillimeter Array (ALMA) 
with its much improved spatial resolution and sensitivity 
can only readily image molecules at '^5 AU or larger in the 
nearest disks. Moreover, ALMA cannot detect molecules 
without a permanent dipole moment such as C2H2 and 
CH4, which are among the most abundant species in the 
inner disk. Results so far show that there is no clear corre- 
lation between the chemistry in the inner and outer parts 
of the disk (Oberg et al. 2011). The chemistry in the inner 
regions seems to be sensitive to different shapes of radi- 
ation fields and the accretion luminosities (|Pascucci et al.l 
I2OO9I : IPontoppidan et al.ll2010D . but these quantities do not 
seem to have an impact on the chemical composition of the 
colder gas further out in the disk. 

Searches for more complex molecules in the inner 
few AU must therefore rely on infrared techniques. 
However, the strong mid-infrared continuum implies very 
low line/continuum ratios for emission lines, even at high 
spectral resolution. Indeed, the recent VLT-CRIRES (i? — 
A/AA = 10^) searches in the 3 /Ltm atmospheric win- 
dow show that lines of molecules other than CO have 
line/continuum ratios of typically only a few %, and that 
even relatively simple species like CH 4 are not detected at 
the ~1% level (|Mandell et al.ll2012[ ). On the other hand, 
absorption lines offer a much better chance of detecting 
minor species for a variety of reasons. First, absorption oc- 
curs from the ground vibrational level where the bulk of 
the population resides, so that the signal is much less sen- 
sitive to temperature. Another advantage is that absorp- 
tion lines are relative in strength to the continuum whereas 
emission lines are absolute. So the strength of the absorp- 
tion lines relative to the continuum will stay the same in 
sources which have a stronger continuum whereas the emis- 
sion lines will be dominated by the continuum. Both these 
advantages imply that absorption lines are easier to detect 
for less abundant molecular species than emission lines. 

Detection of absorption lines requires, however, a special 
orientation of the disk close to edge-on, so that the line of 
sight to the continuum passes through the inner disk. Only 
a few disks have so far been found with such a favorable ge- 
ometry: tlia^jJOundJ)ph-IRS^^ ct al. 2006), GV 



Tau N ( Gibb et al. 2007: jDoppmann et al.l 12008 1 and DG 
Tau B (|Kruger et al. 20111) . In all cases, the mid-infrared 
absorption bands of HCN, C2H2 and CO2 have depths of 5- 
15%, even at the low spectral resolution R « 600 of Spitzer. 
For high S/N > 100 spectra, detection of absorption fea- 
tures of order 1% should be feasible, providing a dynamic 
range of up to an order of magnitude in abundances to 
search for other molecules. With increased spectral resolu- 
tion and sensitivity offered by future mid-infrared instru- 
ments such JWST-MIRI {R «3000), and SOFIA, SPICA 
and ELTs {R > 50,000), another order of dynamic range 
will be opened up. 

A large variety of increasingly sophisticated 
physico-chemi cal models o f the in ner regions o f disks 
exist (e.g., IWillacv et cH 119981: lAikawa et all Il999t 



Markwick et al. 


2002: Nomura et al.l 2007: Agiindez et al. 


20081: Gorti & Hollcnbach' '2008': 'Glassgold et al.' '20091: 
Nomura et al. 2009: WiUacv & Woods 2009: Woitke et al. 


2009; 


Kamo et 


a 


J I2OIO; Walsh et all 


20101: lAresu et al. 


2011; 


Gorti et al.l 


201h: .Heinzeller et al. 


2011l:lNaiita et al. 


2011; 


Vasvunin et al.l2011; Walsh et al. 


20121). The models 



differ in their treatments of radiation fields (UV and/or 
X-rays), the gas heating and resulting disk structure, 
dynamical processes such as accretion flows and disk 
winds, grain properties and chemical networks (e.g., grain 
opacities, treatment of gas-grain chemistry including H2 
formation at high temperature). Some models consider 
only the simplest molecules in the chemistry, others have 
a large chemical network but publish primarily results for 
specie s that can be o bserved at millimeter wavelengths. 
Only iMarkwick et al.l ([^02) list the most abundant 
species, including complex molecules that do not have 
a dipole moment, at 1, 5 and 10 AU as obtained from 
vertically integrated column densities. Since infrared 
observations probe only part of the disk down to where 
the continuum becomes optically thick, these models may 
not be representative of the surface layers. Column den- 
sities appropriate for comparison with i nfrar ed data have 
been p resented by lAgundez et al.l (|2008[ ) and iNaiita et al.l 
(2011) but do not provide data for more complex molecules. 
Woods & WiUacy (2007) and Kress et al. (2010) consider 
PAH processing in the inner disk and study its impact on 
the abundances of related species like benzene and C2H2. 
Note that PAHs are generally not detected in disks around 
T Tauri stars, including the two disks studied here, at 
levels a f actor of 10-100 lower than found in the interstellar 
medium (jGeers et al.ll2006l:IOliveira et al.ll2010D . 

Observations of molecules in comets provide another 
interesting data set for comparison with protoplanetary 
disks. Solar system comets were likely formed at distances 
of about 5-30 AU in the protosolar nebula. Many volatile 
molecules are now routinely observed in cometary atmo- 
spheres at infrared and millimeter wavelengths, including 
spe cies as complex as C2H6, CH3OH, and even (CH20 H)2 
(see lMumma fc Charnlev1l201lHBockelee-Morvanl201ll for 
reviews). It is still debated whether the abundances mea- 
sured in comets directly reflect those found in the dense 
envelopes around protostars or whether they result from 
processing and mixing material from the inner and outer 
disk into the comet-forming zone. Putting constraints on 
the inner disk abundances of these molecules will be im- 
portant to probe the evolution of material from the natal 
protosolar nebula to the formation of icy bodies. 

In this study we use the existing high S/N Spitzer spec- 
tra of IRS 46 and GV Tau to put, for the first time, upper 
limits on various molecules in the inner disk: HNC, CH3, 
C2H4, C2H6, C3H4, HC3N, CgHg, NH3, C4H2, CH3CN, 
H2S and SO2. These molecules were selected to include 
most of the top 15 hi ghest vertical column de nsity molecules 
at 1 and 5 AU by |Markwick et alj (|2002l) . The selected 
species can in principle directly test the predictions of 
models of inner disk chemistry. The list also contains sev- 
eral molecules observed in cometary atmospheres and two 
more molecules with a permanent dipole moment, HNC and 
HC3N, which together with HCN can be observed at both 
infrared and millimeter wavelengths and can thus be used 
to connect the inner and outer disk chemistries through 
ALMA imaging. 
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The mid- infrared spectra of IRS 46 and GVTau con- 
tain detections of C2H2, CO2 and HCN which are analyzed 
here in terms of column densities and abundances, follow- 
ing the same strategy as for IRS 46 in iLahuis et al.l (|2007l) . 
The earlier detections of HCN and C2H2 toward GVTau 
were perfor med in the 3 UTa window with higher spectral 
resolution (jCibb et alJ I2007D but no Spitzer results have 
yet been presented on this source. We have reduced the 
GV Tau spectrum and rereduced the IRS 46 spectrum with 
the latest IRS pipeline in order to allow for a consistent 
comparison between the two sources. For the same reason, 
we have also reanalyzed the IRS 46 spectrum using the same 
analysis software. 

A description of the observations and the reduction of 
the data for IRS 46 and GV Tau is presented in Section [2] 
together with some information about these two protoplan- 
etary disks. Section[3TT]presents an overview of the observed 
and modeled spectra, whereas 13.21 uses a synthetic local 
thermodynamic equilibrium (LTE) model to estimate col- 
umn densities and excitation temperatures for HCN, C2H2 
and CO2 toward GV Tau and compares the results with 
those for IRS 46. Section [5751 presents the upper limits for 
the various molecules toward the two sources and Section 
13.41 shows how future instruments can provide more strin- 
gent limits. Discussion and comparison to chemical models 
is performed in Section |4l A summary of the main conclu- 
sions is found in Section [5l 



2. Observations 

2.1. IRS 46 and GVTau 

The observations of IRS 46 and GV Tau were made using 
Spitzer-mS in both the Short-High (SH; 9.9 - 19.6 /im) 
and Long-High (LH; 18.7 — 37.2 /im) modes with a spectral 
resolving power of i? = X/8X = 600. 

Oph-IRS46 was observed at a = 16''27™29^4 and 
5 = -24°39'16".3 (J2000), located in the Ophiuchus molec- 
ular cloud at a distance of around 120 pc ILoinard et ai] 
120081) ■ IRS 46 was initially observed in 2004 as part of the 
Cores to Disks Spitzer legacy program ()Evans et al.ll2003[ ) 
and in 2008 and 2009 at multiple epochs to search for vari- 
ability. Its mid-infrared spectral energy distribution rises 
stron gly with wavelengt h, as expected for a near edge-on 
disk (iCrapsi et al.l[2008h . Strong HCN, C2H2 and CO2 ab- 
sorption has been detected with Spitzer and attributed to 
arise from warm gas in the surface layers of the inner few 
AU of the disk, seen in absorption against the continuum 
produced by the hot inner rim on the near and far side 
of the star (|Lahuis et al.ll2006[ ). iLahuis et ahl (|201lL and 
in prep.) show that hot water emission lines are also de- 
tected. More interestingly, strong variation in the depth of 
the molecular absorption bands as well as in the strength 
of the water emission lines and the mid-IR contiuum is ob- 
served on timescales of a few years. The data used here are 
the original observations obtained on Augus t 29, 2004 as 
iart o f AOR# 0009829888 and published bv ILahuis et~aLl 
200^. These show the deepest molecular absorptions, thus 
providing the best upper limits of column densities of other 
species relative to the observed C2H2, HCN and CO2 col- 
umn densities. 

GV Tau is a T Tauri star that is partly embedded in 
the L1524 molecular cloud. Its observations were positioned 
at a = 4''29"25^8 and 5 = -h24°33'00"5 (J2000). It 



has an infrared companion about 1.2" to the north. The 
companion is named GV Tau N and the primary optical 
source is called GV Tau S. Gibb et al. (2007) detected HCN 
and C2H2 toward GV Tau N using Keck-NIRSPEC at L- 
band, however, no such detections were made toward GV 
Tau S. GV Tau has subsequently been observed using the 
IRS SH mode at multiple epochs with Spitzer in a G04 
program (PI, F.Lahuis; program ID 50532). For the SH 
part of the spectrum the G04 data (AOR # 0022351616, 
0028247808, 0028247552 and 0031618304) were used. For 
the LH part, data from the Spitzer GTO program ob- 
served on 2 March 2004 as part of AOR # 0003531008 
were adopted. Note that the Spiizer-IRS aperture does 
not resolve the GV Tau binary, in contrast to NIRSPEC. 
The Spitzer spectra therefore combine emission and ab- 
sorption of GV Tau N and GV Tau S. Both GV Tau N 
and GV Tau S are variable at 2 /im (iLeinert et al.l 120011 : 
iKoresko et al.lll999l ). This variability has been attributed to 
variable accretion mechanisms for GV Tau N and variation 
in the extinction due to inhomoge neities in the circumstel- 
lar material for GV Tau S (Leine rt et al.ll200l] ). However, 
the multi-epoch Spitzer data do not show significant mid- 
infrared variation on timescales of a few months up to a 
few years, in contrast with the near-IR ground-based re- 
sults and with IRS 46. The mid-IR continuum of GV Tau 
N is about an order of magnitude brigh ter than the con- 
tinuum of GV Tau S betw een 8 — 13 ^m (|Przvgoddall2004t 
iRoccatadiata et aL[ _ 2011). S i nce n o abs orption was seen 



in GV Tau S in 



Gibb et all ()2007l ) and iDoppmann et al.l 



it is assumed that the majority of the absorption 
arises toward GV Tau N. However, the continuum emis- 
sion from the southern source captured in the Spitzer-UMA 
aperture slightly reduces the total optical depth of the ab- 
sorption lines in the spectrum. To put an upper limit on 
the added uncertainty caused by the additional continuum 
emission from GV Tau S, it is assumed that the mid-IR con- 
tinuum of GV Tau N is ^ 10 times stronger than that of GV 
Tau S, resulting in an additional uncertainty of ~ 1 — 2 % 
for features that are ^ 5 — 20 % deep. Since the effect is mi- 
nor, no correction is made for the column densities derived 
here. 



2.2. Data reduction 

The data reduction started with the BCD (Basic Calibrated 
Data) images from the Spitzer archive processed through 
S18 pipehne. The BCD images were then pr ocessed using 
th e C o res to Disks (c2d) a nalysis pipeline (jLahuis et al.l 
l2006bl : iKessler-Silacci et al.l [2006,)- The mam processing 
steps are background correction, bad-pixel removal, spec- 
tral extraction, defringing, order matching and spectral av- 
eraging. Two extraction methods were used; 1) a full aper- 
ture extraction from the BCD im ages and 2) an opt imal 
extraction using an analytical psf (jLahuis et al.ll2007[ ) de- 
fined using a set of high S/N calibration stars. For both 
extractions a relative spectral response function (RSRF) 
calibration is applied with ^ Dra as the main reference star 
using MARCS mod els taken from the Spitzer science center 
(|Decin et al.ll2004D . 

For all observations the extraction method giving the 
best S/N was used to produce the final spectra. The two 
(partly) independent extraction methods further allow to 
better discriminate between artifacts and true science fea- 
tures. 
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Fig. 1. Spectra of the protoplanetary disks around IRS 46 and GV Tau taken with the Spitzer-IRS. 



3. Results 

3.1. Spectra 

Figure |T] shows the spectra of IRS 46 and GV Tau over the 
10 — 37 /im region. Both spectra differ by more than an 
order of magnitude in the mid-IR continuum level, but they 
show a striking resemblance in spectral shape and spectral 
features. The absorption bands of gaseous C2II2 1^5, HCN 
1^2 and CO2 1^2 can be clearly seen at 13.7, 14.0 and 15.0 fim, 
together with the solid CO2 feature at 15 — 16 ^m. To search 
for other molecules, a local continuum has been fitted to 
the broad spectral features and divided out. The S/N on 
the continuum is typically 100 or better. No other obvious 
absorption features are detected at the few% absorption 
level. The model spectra with derived column densities and 
upper limits are described below. 



3.2. C2H2, HCN and CO2 

To extract quantitative information from the spectra, a 
simple local thermodynamical equilibrium (LTE) absorp- 
tion slab model has been used to fit the data. The free pa- 
rameters in the model are the excitation temperature, the 
integrated column density along the line of sight and the 
intrinsic line width, characterized by the Doppler 6-value. 
The excitation temperature sets the level populations of 
the molecule using the Boltzmann distribution. The lack 
of collisional rate coefficients for many of the species con- 
sidered here prevents non-LTE analyses. The model spec- 
trum is convolved with the spectral resolution of the in- 



strument and resampled to the observed spectra. More de- 
tails about the model and the molecular parameters and 
data that ar e used for the three de t ected n iolecules can 
be found in Lahuis &: van DishoecS (|2000[ ): iLahuis et al.l 
(|2007D : lBoonman et al.l (|2003D and in Table 1. 

Figure ID presents a blow-up of the 13 — 15 /im range 
of the GV Tau and IRS 46 spectra with the continuum di- 
vided out. Included are the best- fitting model spectra. The 
figure clearly shows that the P- and i?-branch lines are 
difficult to detect at the Spitzer-IRS spectral resolution 
of i? = 600, however the Q-branches of C2II2, HCN and 
CO2 are easily seen. In addition the Q-branch changes its 
form and depth with excitation temperature and column 
density. In particular, the depth of the Q-branch decreases 
with increased excitation temperature for the same total 
column density and broadens to the blue side due to an 
increase in the population of the higher rotational levels. 
A higher column density on the other hand increases the 
central depth of the Q-branch si nce more molecules ab- 
sorb photons. As can be seen in iLahuis fc van Dishoecld 
(|200Q) the Q-branch is sensitive to the adopted Doppler 
6-value with the magnitude of the effect depending on 
the temperature and column density. Spectrally resolved 
data obtained wit h Keck-NIRSPEC (Salyk ct al. 201lB 
ILahuis et al] I2006D and within our VLT-CRI RES survey 
(jPontoppidan et all 1201 it iBrown et al.l lsubm.h show that 
the HCN and CO lines have 6 «i 12 km s~^. In our anal- 
ysi s we therefore a dopt a Doppler 6- value of 10 km s~^. 
Lahuis et al.l (|2006t ) used 6=5 km s~^ for IRS 46 which gives 
slightly increased temperatures and reduced column densi- 
ties by ~ 5 - 10%. 
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Molecule 


Formula 


Band 


^1 

[cm J 


c a 
Jlit 

r + -1 -21 
l&tm cm 1 


Q 6 

IcLtni cm 1 


Source^ 


Acetylene 


C2H2 


i/5 CH bending 


729.1 


630 


816 


H08 


Carbon Dioxide 


CO2 


V2 bending 


667.4 


200 


249 


H08 


Hydrogen Cyanide 


HCN 


1/2 bending 


713.5 


257 


286 


H08 


Hydrogen Isocyanide 


HNC 


1/2 bending 


464.2 


1570 


798 


G09 


Methyl Radical 


CH3 


V2 out-of-plane bending 


606.5 


611 


616 


FPH 


Ethylene 


C2H4 


V7 CH2 waggling 


949.2 


324 


320 


G09 


Ammonia 


IN ±13 


U2 symmetric bendng 


yoU.U 


ODO 




xlUo 


Sulphur Dioxide 


SO2 


1/2 bending 


517.6 


113 


97 


G09 


Ethane 


C2H6 


CHS rocking 


822.0 


36 


29 


H08 


Diacetylene / Butadiyne 


C4H2 


vg, CH bending 


627.9 


437 


229 


G09 


Benzene 


CgHs 


CH bending 


673.5 


250 


212 


G09 


Propyne / Methyl Acetylene 


C3H4 


vq CH bending 


638.6 


360 


201 


G09 


Cyanoacetylene / Propynenitrile 


HC3N 




663.4 


278 


94 


G09 


Methyl Cyanide / Acetonitrile 


CH3CN 


Vi CC stretch 


920.3 


6 


3 


G09 



" Central wavenumber of band and band strengths from Constants for molecules of astropysical interest m the gas phase by J. 
Crovisier 

http: / / wwwusr.obspm.fr/~crovisie /basemole / 

' Band strength of the simulated spectra. S — Nl x ^t{v)Su /n, (see App. A.l lHelmichlll99^ . Calculations were performed for 
T = 298 K, 6 = 20kms~^ and n = 1 ■ 10^^ cm~^ at a resolution of 3000 to keep all bands far from saturation 

" H08: HITRAN 2008 (|Rothman et al.ll2009l '). GOO; GEISA 2009 (|Jacauinet-Husson ct al...2011, ') and FPH: iHelmicbl ((l99l l 




13.5 14.0 14.5 15.0 13.5 14.0 14.5 15.0 

Wavelength [/i,m] 

Fig. 2. Continuum normalized spectra of GVTau and IRS 46. Plotted in black are the observed spectra and overplotted 
in red the best-fit synthetic spectra to the absorption bands of C2H2 1^5=1-0, HCN 1/2=^-^ and CO2 i^2=l-0. See Table[5] 
for best fit model parameters. 



Table 2. Results from molecular fits to GV Tau and IRS 46 
absorption features 



Source 


Temperature [K] 
IRS 46 GV Tau 


Column density [lO^*^ cm"^] 
IRS 46 GV Tau 


C2H2 
HCN 
CO2 


490 ± 11 

420 ± i 
250 ± 11 


720 ± |[j 
440 ± II! 
250 ± 11 


2.1 ± 0.4 
3.7 ± 0.8 
8.4 ± 1.1 


1.4 ± 0.3 
1.8 ± 0.4 
5.1 ± 0.7 



A grid of synthetic spectra of C2II2, HCN and CO2 was 
made for a range of column densities and temperatures and 
fitted to the data obtained for IRS 46 and GV Tau. The best 
fit as presented in Fig.[2]was determined by finding the min- 



imum difference between data and model as measured by 
the values. The derived column densities and excitation 
temperatures for the different molecules are summarized in 
T able [21 For IRS 46, the values are consistent with those 
of iLahuis et al.l (j2006ll taking into account that a higher 
6— value was adopted and that the final spectrum presented 
in this paper has a slightly higher S/N ratio compared to 
ILahuis et all (|2006D . It is seen that the temperatures of the 
different molecules and their column density ratios are com- 
parable between GV Tau and IRS 46. This supports the 
hypotheses that both sources are inclined disks with sim- 
ilar characteristics. In both sources CO2 has the highest 
column density but the lowest temperature, whereas the 
HCN/ C2H2 ratio i s sligh tly above unity. In their 2-5 /J,m 
study, iGibb et al] (|2007l) however find significantly lower 
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temperatures for C2H2 (170 ± 20 K) and HCN (115 ± 
20 K) compared to our estimated temperatures of about 
400 to 700 K. However in lat e r Kec k-NI RSPEC L-band ob- 
servat ions. iDopomann et aTl (|2008[ ) and lGibb k, Troutmanl 
()2011|) detect lines out to much higher J values, indicat- 
ing warmer gas around 500 K. Our mid- infrared results are 
therefore not inconsistent with the near-IR data. 

The high resolution near-IR data for IRS 46 and GV 
Tau N show that the spectral lines are shifted in velocity. 
For IRS46, CO and HCN ar e blueshifted by ab out 24 km s'^ 
with respect to the cloud (jLahuis et al.ll2006h whereas for 
GV Tau the HCN lines are redshifted by aboii t 13 km s~^ 
compared with the star (|Doppmann et all2008D . This could 
indicate that the observed HCN, C2H2 and CO2 absorption 
originates in a disk wind or infalling envelope rather than 
the disk itself. However, the high densities needed to ex- 
cite these higher J-transitions as well as the constraints 
on the size of the high abundance region (<11 AU) im- 
ply that the absorption lines have an origin in outflow- 
ing or infalling gas that must be very closely related to 
the disk itself with a chem istry similar to that of the disk. 
This i s further discussed i n Lahuis et al.l (I2006D : iGibb et all 
(I2007D: iDoppmann et all (I2008D: iKruger et al.l (|201lh and 
iMandell et al.l (|2012l ). iFuente et al.l (|2012D recentlv imaged 
the warm HCN associated with the disk of GV Tau N at 
millimeter wavelengths and found an emitting radius of less 
than 12 AU. 



3.3. Other molecules 

3.3.1. Overview and molecular data 

Table |T] summarizes the molecular data (vibrational mode, 
line positions, band strengths) used for all molecules for 
which searches have been made toward IRS 46 and GV 
Tau in the 10-30 /Ltm wavelength range, together with the 
main references from which they have been extracted. Only 
the intrinsically strongest bands of each molecule have been 
chosen; weaker bands are ignored. Note that isotopologues 
and vibrationally excited states or 'hot bands' are not in- 
cluded in the data sets, except in the fitting of the observed 
spectra of HCN, CO2 and C2H2. Hot bands are expected 
to be suppressed in full non-LTE calculations, where the 
excitation of the higher vibrational levels is subthermal at 
densities below « IQ^^ cm~^. Synthetic spectra are gen- 
erated foll owing the procedures a s desc ribed in iHelmichI 
(|1996[) and lLahuis fc van Dishoeckl (|2000f ). 

Our spectra use molecular data from various databases 
listed in Table |TJ To assess their reliability, we have com- 
puted the integrated absorption band strengths of the sim- 
ulated spectra and compared them with independent band 
strengths tabulated in the literature (columns 5 and 6 of 
Table [T]) . This table shows that the band strengths of all 
species agree to within factor of ~ 3 and in most cases 
much better. For the purpose of this paper this accuracy 
is sufficient and differences between databases and other 
literature values are not pursued here. It should be noted, 
however, that there are (sometimes significant) differences 
between line lists and not all line lists are complete. 

Figure |3I presents an overview of the simulated LTE 
spectra for all molecules considered here at T = 500 K, 
R — 50,000 and a column density of 110^^ cm~^. A vari- 
ety of absorption patterns is seen, depending on the charac- 
teristics and symmetries of the individual bands involved. 



Together, they span most of the wavelength range observed 
by Spitzer. 



3.3.2. Upper limits from Spitzer data 

For all species without detections upper limits were derived 
from the local normalized spectrum based on the rms noise 
over a region of 10 resolution elements. The definition of the 
continuum determines the derived rms noise and therefore 
the upper limits. To be conservative, the local continuum 
was defined by a straight line fit to clean regions either 
side of the spectral feature. The upper limit on the column 
density was taken to be the value for which the synthetic 
spectrum has a feature depth of 3cr. Examples for the case 
of GV Tau are presented in Fig. 01 For some species, e.g. 
CH3CN, the model lies below the continuum in the zoomed- 
in region presented in Fig. |4] due to blending of weak P- and 
i?-branch lines at the Spitzer resolution. Over the wider 
wavelength range on which the continuum and model are 
defined both reach unity. 

The depth of a molecular feature depends on the column 
density, the temperature and the adopted Doppler parame- 
ter b. In general, the features become stronger with increas- 
ing column density and decreasing temperature. For high 
column densities, the transitions become optically thick and 
the corresponding absorption lines saturate. The absorption 
depth then no longer increases linearly with column depth, 
with the saturation being stronger for smaller b values. 

Fig. [5] shows how the relative intensity of the spectral 
features change with temperature and column density for 
each molecule at the Spitzer resolution and b = 5kms~^ 
(we adopt a moderately low value for b to illustrate the 
effects of saturation). The effect of saturation is clearly 
seen in some of the spectra, e.g. CH3 at 200 K, where the 
curves deviate from linear at higher column densities. For 
reference. Figures IB. II - IB. 141 in the Appendix present an 
overview of the simulated LTE spectra a.t R = 600 for all 
molecules considered here at additional temperatures of 200 
and 1000 K. Not all molecules follow the expected trend of 
depth versus column density and temperature. This can be 
explained when looking in detail at the different spectra. 
For example, HNC follows the expected trend of decreasing 
depth with increasing temperature fFig. lB.4l) . but not CgHe 
(Fig. IB. lip . The latter behavior is due to the low spectral 
resolving power R = 600 which does not resolve the intrin- 
sically narrow Q— branch of this heavy molecule. At higher 
resolving power, however, the strength of Q-branch does in 
fact decrease with higher temperature as expected. 

Fig.|5]includes the 3ct absorption depth limits for IRS 46 
as dotted lines and the corresponding upper limits on the 
column densities for 200, 500 and 1000 K as dashed lines. 
Saturation is not significant for any of the molecules so 
we can adopt these upper limits for IRS 46 and GV Tau to 
compare with C2H2 and HCN derived with b — lOkms^^. 
The derived 3ct upper limits on the column densities are 
presented in Table [3] and 21 

Tables O and [Bl present the upper limits on the column 
densities relative to C2H2 and HCN, respectively, for both 
IRS 46 and GV Tau. Abundance ratios relative to C2H2 and 
HCN are typically of order unity, except for CH3CN and 
C2H6 which have particularly low band strengths (see Table 
[T]). The most stringent ratios of < 0.2 — 0.5 are obtained for 
C4H2 and CgHg. The ratios are also graphically displayed 
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Fig. 3. Synthetic spectra of the different molecules for Tcx^SOO K and for a resolving power of i?=50,000. The column 
density of each molecule is set to be 110^^ cm^^ which is the same as in Fig. IB.1IIB.13I 



in Figures [6] and [7] where they are compared with model 
results. These results will be further discussed in fJH 



3.4. High resolution spectra 



Figures IB.li IB. 21 and IB. 31 in the Appendix present C2H2, 
HCN and CO2 absorption spectra at higher spectral res- 
olution for R ss 3000, as appropriate for the JWST- 
MIRI instrument, and at i? ~ 50,000, as typical for fu- 
ture mid-infrared spectrometers on an Extremely Large 
Telescope (ELT). The latter spectrum is also character- 
istic (within a factor of 2) of the s pectral re solution of 
R = 100, 000 of EXES on SOFIA (Ri chter et al.l 20061 or a 
pote ntial high resolution spectrome ter on the SPICA mis- 
sion (jGoicoechea fc Nakagawal[201lt) . As expected, the cen- 



tral Q branch becomes deeper with higher spectral resolu- 
tion and the P— and i?-branches become readily detectable, 
allowing a more accurate model fit to the data. The inferred 
column density from such data should not change beyond 
the error bars derived from the low resolution data, how- 



Figures IB. 41 - IB. 121 present the higher resolution spec- 
tra of all other molecules at i?=3000 and 50,000 at Tox — 
200, 500 and 1000 K using 6 = 5 km s~^ and a column 
density of 110^^ cm~^. The improved detectability of the 
molecules at higher spectral resolution is obvious. To illus- 
trate the improvements compared with the current Spitzer- 
IRS observations, the upper limits for the high resolution 
cases are included in Table |3] for IRS 46 assuming the same 
S/N values. For GV Tau, the limits at higher resolution 
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Fig. 4. Blow-ups of synthetic spectra (in red) for different molecules at a 3cr maximum optical depth compared with the 
observed spectrum of GV Tau (in black) . 



scale similarly. The column density limits are lower by fac- 
tors of 2-10. They do not decrease linearly with increasing 
resolving power, however, because the strong Q— branches 
used to set the limits are blends of many lines at low resolu- 
tion which become separated at higher resolving power. In 
some cases (C2H6, CH3CN), the gain is very small because 
the molecule absorbs less than 0.1% of the continuum for 
the adopted column density, which remains undetectable at 
S'/iV=100 even at high spectral resolution. 

Tables [S] and IH] include the column density ratios with 
respect to C2H2 and HCN at higher resolving power. Since 
the C2H2 and HCN column densities remain the same, the 
abundance ratio limits are now up to an order of magnitude 



lower, thus bringing the limits in a more interesting regime 
where they provide more stringent tests of chemical models. 
To push the abundance ratios even lower, a higher S/N 
than 100 on the continuum is needed. 

Note that these limits do not take the transmission 
of the Earth's atmosphere into account but assume that 
the strong Q-branches can be observed unobscured. For 
ground-based instruments and to a lesser extent SOFIA, 
this is often not the case and detectability depends both 
on the transmission and on the radial ve locity sh ifts of the 
sources with respect to atmospheric lines (jLacv et al. 198^. 
For JWST, the better stability of a space-based platform as 
well as the large mirror and the absence of telluric absorp- 
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Fig. 5. Variation of the maximum depth of the spectral features as a function of column density at excitation temperatures 
of 200 K (red), 500 K (blue) and 1000 K (black). The 3cr observational limit for IRS 46 is marked with a black dotted line 
and the colored dashed lines show the corresponding upper limits on the column densities for the three temperatures. 



tion bring significantly higher S/N ratios than 100 within 
reach once the instrument performance is well characterized 
in orbit. 

4. Discussion 

Tables [5] and [6] and Figures [6] and [7] compare our limits with 
a variety of chemical models. There are two distinct routes 
towards molecular complexity in regions of star- and planet 
formation. First, at elevated temperatures such as found in 
the inner disks, various reactions with activation barriers 
open up. If atomic carbon can be liberated from CO and 



atomic nitrogen from N2, high abundances of CH4, C2H2 
and HCN can be produced. The second route starts in the 
pre-stellar cores where ices are formed through grain surface 
reactions. At a later stage, these ices can be transported 
into the disk and evaporate so that a chemistry rich in 
hydrogenated molecules can ensue. We review each of these 
classes of models and then discuss our observations in the 
light of these models and in comparison with cometary and 
other data. 
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Fig. 6. Comparison of observed inner disk abundance ratios (wit. C2H2) with chemical models (upper panel) and 
cometary observations (lower panel) for IRS 46 and GVTau. In both panels detections (diamonds) and upper limits 
(arrows) for IR S 46 ( black) and GV Tau (red) are indicated. Upper panel: Abundance ra tios in the di s k mo del by 
iMarkwick et alJ (2002) at 5 AU (green triangle), from the reference disk mo del a t 1 AU bv iNaiita et al.l (j201lt) (blue 
square) and at a 0/C = 1 (blue cross), and from the disk model bv .Agundez et al.l ( 20081) at 1 AU (green s quare ) and 3 
AU (green cross). Lower panel: Observed range of cometary abundance ratios from Mumma fc Charnlev I (|2011[ ) (green 
bar and stars). 



4.1. Warm chemistry 



The warm gas chemistry in the photospheres of disks follows 
a similar chemical scheme as that in other interstellar re- 
gions with high temperature gas such as t h e inner envelopes 



of massive protostars (e.g., Lahuis fc v an Dishocck* 1200' 



Doty et all [20021: iRodgers fc Charnl'evir2003.; Stauber et al.l 
20051) ■ high density photodiss ociation regions (PD Rs) (e.g., 
Sternberg fc Dalgarnol r995fl and shocks (e.s... | Mitchell| 



1984; Pincau dcs Fo rets et al.lll987t IViti et al.lboill ). The 
chemical scheme starts with separating C or C"*" from CO 
and N from N2, see Fig. |8] (adapted from lAgundez et al.l 
[2OO8I) . This can be done either by UV photons, cosmic rays 
or X-rays. A UV-dominated region will produce compara- 
ble amounts of C^ and C, whereas cosmic rays and X-rays 
produce He^ which will react with CO to produce primarily 
C+. 



For carbon-bearing species, the rate limiting steps of 
this scheme are the reactions of C and C+ with II2 which 
have activation barriers Ea of ~12,000 and 4,000 K, re- 
spectively. At temperatures of a few hundred K, the C+ 
channel leads to CII4 and C2H2 while at high tempera- 
tures (>800 K), the C-channel becomes active. Reactions 
of C2 with H2 have activation barriers of ^^1500 K and 
form another route to produce C2H and subsequently C2II2. 
For nitrogen-bearing species, the reactions of NH with H2 
{Ea ~ 7800 K) to form NH3 and of CN with H2 (Ea - 820 
K) to form HCN also require high temperatures. At low 
temperatures < 200 K, the above reactions are closed and 
C+, C and N will be driven back to CO and N2 through 
reactions involving OH. Hence, the higher the temperature 
the more C2H2, CH4 and HCN will be produced. Since the 
temperature decreases with disk radius, the abundances of 
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Fig. 7. Comparison of observed inner disk abundance ratios (wrt. HCN) with chemical models (upper panel) and cometary 
observations (lower panel) for IRS 46 and GV Tau. Use of symbols and references identical to Figure IHl 



C2H2, CH4 and HCN show a st rong radial dependenc e with 
a steep decrease beyond 1 AU (jAgundez et al.l[2008h . 

Similar arguments apply to the more complex hydrocar- 
bons studied here. In fact, high temperature, high density 
chemistry starting with high abundances of C2H2 and HCN 
resembles the chemistry of the atmospheres of carbon-rich 
evolved stars, which has been studied for decades (e.g., 
ICherchneff fc GlassgoldlfTool iMillar fc Herbsd [1991 . For 
example, C2H2 reacts with CH3 produce C3H4, with C2H 
to C4H2, and with NH3 to CH3CN. If C2H is sufficiently 
abundant (due to photodissociation of C2H2, for example), 
subsequent reactions may lead to large unsaturated car- 
bon chains. More saturated molecules such as C2H4 and 
C2Hg need CH4 as their starting point. Finally, nitrogen- 
containing species are produced by reactions of HCN or N 
with hydrocarbons. For example, C2H2 and C2H react with 
HCN to form HC3N. It is now well established that many 
of these neutral-neutral reactions actually have substantial 
rates even at low temperatures (<100 K) ((Smith 2011.'). 
Thus, the bottleneck in producing the more complex or- 
ganic molecules is the formation of the chemical precursors 



C2H2, HCN and CH4. The more complex molecules are thus 
also expected to show a strong dependence on disk radius. 

Benzene is the most complex molecule probed here and 
is a precursor for building larger PAHs. One route to form 
benzene is through the reaction of C3H4 and its ion, found 
to be important in the in ner disk chemistry in some models 
(jWoods fc Willacvl[2007[) . Alternatively, the destruction of 
PAHs by high temperature gas-phase reactions inside the 
'soot line' can lead to the production of small hydrocarbons, 
including C2H2 and perhaps CgHg ([Kress e t al. 2010). 

CO2 is the product of the reaction of OH with CO. 
Models show that all atomic O is driven into H2O at tem- 
peratures above ~230 K if th ere arc n o rapid H9O de - 
struction routes (e.g., iDraine e t al. 1981; ICharnlevlll997( ). 
The formation of OH also needs elevated temperatures, 
but peaks in the 100-200 K range because of the rapid 
reactions of OH with H2 to form H2O at higher temper- 
atures. In warm gas, the CO2 abundance follows that of 
OH and also peaks at 100-200 K, with a strong decrease 
toward higher temperatures. Thus, in contrast with C2H2 
and HCN, the abundance of CO2 increases somewhat be- 
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Table 3. Inferred upper limits of column densities [10^^ cm^^] toward IRS 46 at different excitation temperatures and 
spectral resolving powers. 







C2H4 


C2H6 


CeHe 


CH3 


HNC 


C3H4 




=600, 200 K 


<5.1 


<69 


<0.7 


<0.6 


<0.6 


<2.0 




500 K 


<3.2 


<110 


<0.5 


<r^ 6 




<6.2 




1000 K 


<4.3 


<170 


<0.5 


<2.6 


<2.2 


<13 


R= 




<1.8 


<50 


<0.2 


<0.4 


<0.9 


<3.1 


R= 


^50,000, 500l<°- 


<0.6 


<29 


<0.09 


<0.06 


<0.09 


<1.8 


3(7 


1%]'' 


3 


3 


2 


2 


5 


2 






HC3N 


C4H2 


CH3CN 


SO2 


NH3 






=600, 200 K 


<1.1 


<0.7 


<660 


<19 


<2.0 






500 K 


<1.2 


<0.7 


<780 


<25 


<3.9 






1000 K 


<1.4 


<0.9 


<180 


<46 


<7.7 




R 


- 3000, 500 K 


<0.4 


<0.3 


<310 


<8.5 


<14 




R 


- 50,000, 500 K 


<0.3 


<0.2 


<46 


<2.6 


<0.3 




3a 




2 


3 


3 


2 


3 





" Upper limits at spectral resolving powers of _R = 3000 and 
50,000 at r = 500 K adopting S/N = 100, similar to the current 
Spitzer-IRS data at lower spectral resolution. 

' The 3(7 limit in % absorption at the location of the molecular 
band as measured in the Spitzer-IKS data. 



Table 4. Inferred upper limits of column densities [10^^ cm ^] toward GV Tau at different excitation temperatures at 
a resolving power of i?=600. 





C2H4 


C2H6 




CH3 


HNC 


C3H4 


200 K 


<1.2 


<23 


<0.3 


<0.4 


<0.6 


<2.0 


500 K 


<0.9 


<38 


<0.3 


<1.2 


<1.2 


<6.2 


1000 K 


<0.9 


<57 


<0.3 


<1.8 


<2.2 


<13 


3a [%\ ^ 


1 


1 


1 


1.5 


5 


2 




HC3N 


C4H2 


CH3CN 


SO2 


NH3 




200 K 


<1.1 


<0.2 


<230 


<14 


<0.9 




500 K 


<1.2 


<0.2 


<260 


<19 


<1.9 




1000 K 


<1.4 


<0.3 


<50 


<34 


<3.7 




3a [%] 


2 


1 


1 


1.5 


1.5 





" The actual 3cr limit in % absorption at the location of the 
molecular band. 



yond 1 AU in disk models, u ntil the temperature falls below 
100 K (lAgundez et al.ll2008l^ . 

Another parameter that can affect the chemistry is 
the 0/C ratio in the gas. There have been suggestions 
that this ratio changes with disk radius due to migra- 
tion of i cy planetesimals conta ining a large fraction of the 
oxygen (jCiesla fc Cuzzil l2006l ). Hydroc arbon abundances 
are very sensitive to the 0/C ratio (iLanger fc Graedel 
ll989HBettens et al.lll995l: iNaiita et al.ll201lf ) As expected, 
oxygen-bearing species such as H2O and O2 decrease dra- 
matically in abundance when 0/C is less than 1 because 
all volatile oxygen (rather than all carbon) is locked up in 
CO. Pure carbon species such as C2H2 increase in abun- 
dance by orders of magnitude while nitrogen species only 
change appreciably when they contain C (e.g., HCN). Since 
both IRS 46 and GV Tau show large columns of gaseous 
CO2, the lines of sight through the disks probed by the 
mid-infrared data must pass through gas with 0/C> 1. 



Finally, the overall gas /dust ratio c an play a role (e.g., 
INaiita et all201lHBruderer et al.ll2012l) . However, compar- 
ison of the CO column density with the extinction measured 
from the silicate optical depth shows that this ratio is close 
to the interstellar value for both IRS 46 and GV Tau, in- 
dicating no significant grain grow th and settling alon g the 
lines of sight through these disks (jKrueer et al.ll201]] ). 

In SjX]in the Appendix, more detailed descriptions of a 
few specific protoplanetary disk models are given, focussing 
on the simpler species which can be compared with obser- 
vations. 



4.2. Surface chemistry 

The more complex organic molecules such as dimethyl 
ether seen toward protostars are the product of an ac- 
tive gas-grain interaction (se e iTielens fc CharnlevI 119971: 
iHerbst fc van Dishoeckl l2009l for reviews). Recent mod- 
els produce these species on the grains rather than in 
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Table 5. Observed molecular column density ratios relative to C2H2 in disks compared with different chemical models 
and cometary observations. 



Ratio rel to C2H2 


C2H4 


C2H6 




CH3 


HNC 


C3H4 


HC3N 


CH4 


1Kb 4d 


<1.0 


<02 


<0.2 


<U.8 


<U.D 


^ Q 

<3.U 


<U.D 




GV Tau 


<0.6 


<27 


<0.2 


<0.9 


<0.9 


<4.4 


<0.9 


- 


DR Tau 


- 


- 


- 


- 


- 


- 


- 


<2.3^ 


R = 3000^ 


<0.9 




<0.08 


<0.2 


<0.4 


<1.5 


<0.2 




R = 50, 000'' 


<0.3 


<H 


<0.04 


<0.03 


<0.04 


<0.9 


<0.1 




Markwick et al. (2002) 1 AU 








1.0 


1.0 


1.0 


1.0 


1.1 


Markwick et al. (2002) 5 AU 










1.0 


1.0 


1.0 


1.1 


Naiita ct al. (2011) ref 








4.0E-4 








0.09 


Naiita et al. (2011) 0/C = 1.0 








2.5E-6 








3.2E-5 


AEundez et al. (2008) 1 AU 










0.03 






0.01 


Aeundez et al. f2008~) 3 AU 










0.03 






0.2 


Mumma & Charnlev (2011) 




0.2-20 






0.01-0.5 




0.01-0.7 


0.2-15 


Ratio rel to C2H2 


HCN 


C4H2 


CH3CN 


SO2 


NH3 


CO2 







IRS 46 


1.8 


<0.3 


<370 


<11.9 


<1.9 


4.0 


GV Tau 


1.3 


<0.2 


<190 


<13.6 


<1.4 


3.6 


DR Tau 














R = 3000'' 




<0.1 


<150 


<10 


<0.7 




R = 50, OOO*" 




<0.1 


<22 


<1.2 


<0.1 




Markwick et al. (2002) 1 AU 


1.0 


1.0 


1.0 


0.9 


1.0 


1.0 


Markwick ct al. (2002) 5 AU 


1.0 


1.0 


1.0 


0.9 


1.1 


1.1 


Naiita et al. (2011) rcf 


70 






8700 


310 


110 


Naiita et al. (2011) O/C = 1.0 


0.9 






1.1 


0.01 


0.06 


AEundez et al. (2008) 1 AU 


3.0 








2.0E-4 


20 


Aeundez et al. (2008) 3 AU 


1.0 








0.03 


2000 


Mumma & Charnlev (2011) 


0.2-5.0 




0.02-0.4 


0.2-1.0 


0.6-20 


5.0-300 



" From 3 /im emission data (|Mandell et al.ll2012D . 

' The i?=3000 and 50,000 values refer to the Umits on the 
abundance ratios that can be obtained toward IRS 46 if observed 
at higher spectral resolution at a similar S/N. 



hot gas-phase chemistry, either through direct grain sur- 
face chemistry or through mild photolysis of simple ice 
species resulting in radicals which th en react with each 
other to form more complex species ()Garrod et al.l [2008t 
lOberg et al]|2009f) . These types of grain chemistry - either 
direct surface reactions or photolysis - will not result in 
high abundances of C2H2 and HCN because these species 
will be readily hydrogenated by H. On grain surfaces, N is 
quickly converted to NH3, while in the gas phase, N flows 
through NO to N2 with a slight detour to HCN. Grain 
surface reactions lead to very high abundances of NH3 , es- 
pecially in relation to HCN. Observations and comparisons 
with model predictions show that also CH4 formation is 
very efficient on dust grains (|Oberg et al.l 120081) . 

4.3. Comparison of models with observations 

Table [5] and Figures Eland [7] show that warm chemistry can 
explain the observed abundance ratio s of HCN, C2H2 and 
CO2 as well as most limits. In the M arkwick et al.l (|2002[ ) 
models, all abundance ratios are close to unity, but these ra- 
tios refer to the entire disk rather than just the surface lay- 
ers. It is therefore not possible to p roper ly tes t these mod- 
els. T he models of lAgiindez et al.l (2cio8') and ' Naiita et al.l 
(j201l|) provide column densities for just the warm surface 
layers probed by the mid-infrared data, but tabulate only 
a limited number of species see ^|2]for details. 



The inner disk models of lAgundez et al.l ()2008[ ) provide 
the correct range of column densities and abundance ra- 
tios in the surface layers around 1 AU. At 3 AU, the C2H2 
column has dropped dramatically due to the lower tem- 
perature and the model overproduces the CO2/C2H2 ratio 
by a factor >100. The lower excitation temperatures ob- 
served for CO2 in our data compared with C2H2 and HCN 
provide indirect evidence for the radial dependence of the 
chemistry, with CO2 peaking at larger distances (see S 14. II) . 
The pure XDR model results given bv lNaiita etakl ()201lD . 
for the reference case of a disk at 1 AU, agree less well 
with the observations than the PDR model results at 1 AU. 
For example the SO2, NH3 and CO2 model ratios relative 
to C2H2 are several orders of magnitude higher than our 
observed upper limits and det ections. This can be largely 
explained by the fact that the iNaiita et al.l ()2011[ ) models 
do not contain photodissociation reactions in the surface 
layers which limit their abundances. Lowering the 0/C ra- 
tio to 1 instead of 2.5 decreases the difference between the 
observations and the pure X- ray model s but n ow underpro- 
duces CO2. As shown by Wal sh et al.l (|2012[ ) in combined 
UV -I- X-ray models, the inclusion of UV-induced processes 
is generally more important than those due to X-rays. 

The importance of gas-grain interaction for the compo- 
sition of the inner regions of the disk may well be revealed 
by NH3 and CH4 searches. As mentioned in section l4?2l high 
abundances of CH4 and NH3 relative to HCN and C2H2 
could indicate that the chemistry in the inner part of disks 
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Table 6. Observed molecular column density ratios relative to HCN as well as CH3CN/NH3 compared with chemical 
models and with cometary observations. 



Obs/models 


HNC 


HC3N 


CH3CN 


NH3 


CH3CN/NH3 


1Kb 4d 


<U.3 


<U.O 


^ 1 n 
<210 


^1 1 
<1.1 


<200 


GV Tan 


<0.7 


<0.7 


<140 


<1.1 


<140 


AS 205 


- 


- 


- 


<2.7^ 


- 


R = 3000^ 


<0.2 


<0.1 


<84 


<o.r 


- 


R = 50, OOO*" 


<0.03 


<0.08 


<12 


<0.08 




Markwick et al. (2002) at 1 AU 


1.0 


1.0 


1.0 


1.0 


1.0 


Markwick et al. (2002) at 5 AU 


1.0 


1.0 


1.0 


1.1 


0.9 


Naiita ct al (2011) ref 






4.0 






Naiita et aT (2011) 0/C = 1.0 






0.01 






AEundez et al. (2008) 1 AU 


0.01 






6.0e-5 




Aeundez et al. (20081 3 AU 


0.03 






0.03 




Mumma & Charnlev (2011) 


0.01-0.6 


0.01-0.9 


0.02-0.5 


0.6-25 


0.01-0.1 


Garrod et al. (2008~l M 








890 


7.5E-5 



" From 3 /im emission l|Mandell et al.|[20ll ). 

The i?=3000 and 50,000 values refer to the 3(j limits on the 
abundance ratios that can be obtained toward IRS 46 if observed 
at higher spectral resolution at a similar S/N. 



has been reset due to evaporation of ices fror n dust grains. 
Table ini includes the results from model M of iGarrod et al.l 
(|2008l) . The efficiency to form NH3 on dust grains is the 
main reason why NH3 /HCN and CH .3CN/NH3 are so dif- 
ferent in the grain surface models of iGarrod et al.l ()20Q8[ ) 
compared w it h the warr n chemistry models d eveloped by 
iNaiita et all (|201lf ) and lAgundez et al.l (|2008D . Currently 
no observational limit can be set on the CH3CN/NH3 ratio. 
However, our observed limit on NH3/HCN clearly favors a 
low ratio for the abundances of these species, inconsistent 
with a significant contribution from pure ice chemistry. 

The observed upper limits on the other species inves- 
tigated here are in general higher or comparable to what 
the chemical models predict. Thus, a firm conclusion from 
our data is that the molecular abundance ratios cannot 
be higher than what is predicted in the current models. 
However they can be lower and future higher resolution 
observations with for example JWST or ELTs that go an 
order of magnitude (or more) deeper may be able to di- 
rectly test the chemical models (see Table [S]). Specifically, 
deep searches for HNC, HC3N, CeHg, SO2 and NH3 should 
be able to distinguish between some of the models, espe- 
cially if the S/N ratio of the data can be pushed to values 
significantly higher than 100. 



4.4. Comparison with protostars, other disks and comets 

The abundance ratios of the detected molecules — HCN, 
C2H2 and CO2 — are remarkably close (within factors 
of two) to those observed to ward high-mass protostars 
(|Lahuis &: van Dishoeckl |2000[ ). which have been inter- 
preted with high temperature gas-phase chemistry models 
(|Dotv et al.ll2002D . Since our NH3/HCN limits of < 1.1 are 
close to the detected NH3/HCN ratio of ~1.2 toward one 
high- mass protostar (jKnez et al.ll2009[ ). this suggests that 
deeper high resolution searches for NH3 may be fruitful. 

It is also interesting to compare the two observed 
C2H2/HCN abundance ratios of 0.6 and 0.8 for IRS 46 
and GV Tau respectivel y with the one order o f magnitude 
lower values presented in lCarr fc Naiital (|201l[ ) for 5 proto- 



planetary disks in Taurus. Those ratios have been inferred 
from emission lines, whereas ours come from absorption 
data and may probe a different part of the disk (< 1 AU 
vs up to a few AU). This could imply that the C2H2/HCN 
abundance ratios probed through our edge-on disks pre- 
sented ar e more similar t o the hot cores and comets in 
Fig. 12 of |C arr fc Naiital (l201ll) than t o th eir protoplane- 
tary disks. However, 'Sal vk et al.l (l20TTh and iMandeh eTal! 
(2012) present C2H2/HCN abundance ratios of '^l for sev- 
eral protoplanetary disks, illustrating that different anal- 
ysis methods of the emission data can lead to an order of 
magnitude different abundance ratios. The latter values are 
about the same as the abundance ratios found in the hot 
core observations and in our disk absorption data. 



It is interesting to further compare our inferred molec- 
ular ratios to those observed in comets, see Table |5| and 
Figures [5] and [71 The ratios tabulated for comets are pre- 
sented as ranges between the high est and lowest observed 
ratios umma fc Charnlev I [20ll . and are all within our 
observed upper limits and detections. This comparison 
suggests that some of the HCN and C2H2 produced by 
warm chemistry in the inner disk may be incorporated into 
cold comets. However, the presence of complex organics 
such as dimethyl ether in high abundances in comets (see 
iBockelee-MorvanI |2011L for a recent review) suggests that 
the ice chemistry route is also important for their organic 
inventory. The detection of additional species in the inner 
disk such as CH4, HNC and HC3N will further constrain 
their relative importance. 



In summary, our observed ratios suggest that warm 
chemistry models with 0/C>l are most relevant for ex- 
plaining the observed abundance patterns in disks. Future 
higher resolution observations of these molecules with 
JWST-MIRI and other facilities can help in answering the 
question to what extent warm chemistry and surface chem- 
istry contribute to the chemical composition of the gas in 
the planet- and comet-forming zones of disks. 
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Fig. 8. Chemical diagram describing the main reaction routes to form HCN, C2H2, CH4 and NH3 in warm gas in 
protoplanetar y disks or hot cores d epending on the radiation fields, temperature and different molecular abundances, 
adapted from lAgundez et all (|2008f ). At high temperatures of several hundred K the reaction routes (in red) starting 
with C reacting with H2 to CH and NH to NH2 dominate over lower temperature chemical routes (in blue) of C+ going 
to CH+ and NHj^ and NH^. The different activation energy (Ea) barriers for the chemical reactions are represented by 
different types of arrows; the thicker the arrow, the higher the activation energy. 



5. Conclusions 

In this paper, mid-infrared spectra of HCN, C2H2 and CO2 
have been analyzed for two edge-on disks, IRS 46 and GV 
Tau. The high S/N data have also been used to put upper 
limits on the abundances of other molecules predicted to 
be abundant in the inner disk. The main conclusions are: 

• The two disks have similar column densities and similar 
abundance ratios of warm HCN, C2H2 and CO2. The 
first two molecules probe gas with excitation tempera- 



tures T ^ 400 — 700 K, whereas CO2 probes somewhat 
cooler gas. These results are similar to those found to- 
ward massive protostars. 
• The observed abundance ratios for these three molecules 
are consistent with high temperature inner disk chem- 
istry models with 0/C>l in which the abundances of 
HCN and C2H2 rapidly drop with radius beyond 1 AU 
due to the decrease in temperature but that of CO2 
peaks at larger radius. 
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• No other absorption features are detected above 3cr in 
either source, providing upper Umits on a variety of hy- 
drocarbon molecules, NH3 and SO2 that are of order 

unity or less with respect to C2H2 or HCN. 

• The upper limits relative to C2H2 and HCN are ci- 
ther higher or close to values given by high temperature 
chemistry models of protoplanetary disks. The observed 
NH3 /HCN limit is much lower than would be expected 
if the chemistry in disks would have been reset due to 
evaporation of icy mantles on dust grains. 

• Hot chemistry disk models including both UV radiation 
and X-rays produce abundance ratios in better agree- 
ment with our observations than pure X-ray models. 

• The observed abundance ratios in comets arc within the 
same range as our observed ratios or upper limits. The 
composition of comets could therefore be partly build 
up from gas in the inner regions of protoplanetary disks 
mixed outward to the comet-forming zone. 

• Future observations using higher resolution instrumcints 
on JWST, ELT, SPICA and SOFIA will be able to de- 
tect column densities which are an order of magnitude 
or more lower than the upper limits extracted from the 
Spitzer-IRS data, especially if S/N ratios much higher 
than 100 can be obtained. Such data would provide 
much better constraints of the hot gas phase chemical 
models of the inner disk. Edge-on systems such as IRS 
46 and GV Tau remain uniquely suited for this purpose. 
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Appendix A: Discussion of specific disk models 

In this section, we provide more details of the specific disk 
models to which our observational results are compared. 



A.l. X-ray dominated region (XDR) surface layers 



iNaiita et al.l (|2011[ ) have analyzed the chemistry of the in- 
ner (0.25-20 AU) portions of protoplanetary disks exposed 
to X rays. The physical struct ure derives from the th ermo- 
chemical model developed bv iGlassgold et all (|2009[ ). The 
thermal structure is evaluated separately for dust and gas 
and the surface gas temperature significantly exceeds that 
of the dust in the upper layers. In their model, gas and 
dust temperatures are decoupled for column densities less 
than 10^^ cm~^. If only the disk surface is considered (N < 
3- 10^^ cm~^), gas temperatures of 300 K are reached out 
to radial distances of 4 AU. 

The formation of hydrocarbons in this model is trig- 
gered by X-rays which produce He"*" which liberates the 
0+ from CO and N from N2. Acetylene is subsequently 
produced by C and C"*" insertion reactions with small hy- 
drocarbon radicals (Fig. [SJ. Hence, the C2H2 abundance 
is sensitive to the X-ray ionization rate. Radicals such as 
OH are also sensitive to X-rays. At low temperature, OH 
results from the recombination of HsO^ produced by ion- 
molecule reactions and its abundance scales directly with 
the ionization rate. The formation rate of daughter species 
of OH such as NO, SO, SO2 and CO2 is similarly increased 
by X-ray ionization. However, if their destruction is also 
dominated by He^, their abundances are not sensitive to 
the X-ray ionization flux. This is the case fo r CO2 , for ex- 
ample. Note that the models of lNaiita et al.l (|201l[ ) do not 
include UV photodissociation and may thus underestimate 
the amount of OH that could be produced from H2O in the 
upper layers. 

The transformation of N back to N2 is mediated by 
neutral reactions and is accompanied by significant column 
densities of warm nitrogen-bearing molecules such as HCN 
and NH3 in these models. Specifically, HCN is infiuenced 
by X-rays in its formation route both because it liberates 
N from N2 and then N reacts with OH to form NO, with 
NO subsequently reacting with C to form CN. HCN is then 
formed through reactions of CN with H2 . Because of lack of 
photodissociation in the model, HCN is destroyed by He+ 
and hence destruction is also sensitive to the X-ray. In the 
end, the HCN decreases slightly w ith X-ray luminosity in 
the models bv lNaiita et all (|201lh . 



high abundances of C2H2 and CH4 requires the reaction of 
atomic C with H2 to proceed which has a very high acti- 
vation barrier. This reaction only proceeds in very warm 
gas (T > 500 K) forming CH. This is reflected in the radial 
dependence of the C2H2 and CH4 columns, which reach 
values of '^10^^ cm^^ out to radn of ~1 AU, but then drop 
by orders of magnitude in the colder gas. This rapid drop 
with disk radius is also seen in models by other authors. 
In contrast, the CO2 column increases with radius in the 
in ner 1 AU sinc e it fav ors somewhat colder gas. 

IWalsh et al.l (|2012f ) present a combined UV -I- X-ray 
model which shows that a correct treatment of the pho- 
todissociation is generally more important than including 
X-rays. For the important species considered here — C2H2, 
HCN, CO2, CH4 and NH3 — the column densities in the in- 
ner disk do not change measurably when X-rays are added 
to the UV model. 



Appendix B: Auxilliary figures 

This appendix presents simulations of the spectra of all 
molecules considered here at higher spectral resolving 
power of i? = 3000 and R = 50, 000, appropriate for future 
instruments. In addition, spectra at the Spitzer resolving 
power of i? = 600 are included. All spectra are computed 
for Tex = 200, 500 and 1000 K, 6 5 km s'^ and a column 
density of 110^^ cm^^. 



A. 2. Photodissociation region (PDR) surface layers 

Various recent models have analyzed the gas phase 
chemistry of hot inner regions of protoplanetary disks 
including UV radiation for the chemistry and heat- 
ing of the gas (lAeundez et al.l 120081: IWoitke et al.l l2009t 



Willacv fc W oods"2009': 'Vasvunin et all 1201 it IWalsh et alJ 
2012f) ThelAgundcz et al. (2008) study is particularly in- 



structive because they present models for just the photo- 
sphere of the disk, down to H2 column densities of 5 x 10^^ 
cm^^ to which the UV penetrates. In their models, FUV 
photons and cosmic rays produce C, C"*" and N. Atomic N 
is then channeled to HCN through a similar reaction rou- 
tine as that described in Fig. [S] The C"*" leads to a rich 
hydrocarbon radical chemistry. However, the formation of 



Jeanette E. Bast et al.: Exploring organic chemistry in planet-forming zones, Online Material p 2 



>< 1.00 

3 



T3 

!U 0.96 
N 

P 

^ 0.92 
O 

2 1.00 



0.92 



T = 200 K 



R ~ 600 



Tirffffirt 



T = 500 K 



T = 1000 K 



12 



13 



15 



16 



1.00 
0.95 
0.90 
0.85 
1.00 
0.95 
0.90 
0.85 
1.00 
0.95 
0.90 
0.85 



- 


R ~ 3000 






, '■'''■mT|'|'[i|'|ii|i|i|y|i^ 'Tr|'|i|i|i|i|i|i|||'nTiTi'' 



12 13 14 15 16 

Wavelength [yum] 



1.0 



0.6 



0.6 



0.6 



12 



13 



R ~ 50,000 



15 



16 



Fig. B.l. The synthetic spectrum of C2H2 at a column density of l.OTO^^ cm ^, excitation temperatures of 200 (top), 
500 (middle) and 1000 K (bottom), and spectral resolving powers of 600 (left), 3000 (middle) and 50,000 (right). Note 
the different vertical scales for the different spectral resolving powers in this and subsequent figures. 
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Fig. B.2. The synthetic spectrum of CO2 at a column density of 1.010^^ cm ^, excitation temperatures of 200 (top), 
500 (middle) and 1000 K (bottom), and spectral resolving powers of 600 (left), 3000 (middle) and 50,000 (right). 
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Fig. B.3. The synthetic spectrum of HCN at a column density of 1.010^^ cm ^, excitation temperatures of 200 (top), 
500 (middle) and 1000 K (bottom), and spectral resolving powers of 600 (left), 3000 (middle) and 50,000 (right). 
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Fig. B.5. The synthetic spectrum of CH3 at a column density of 1.010^^ cm ^, excitation temperatures of 200 (top), 
500 (middle) and 1000 K (bottom), and spectral resolving powers of 600 (left), 3000 (middle) and 50,000 (right). 
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Fig. B.6. The synthetic spectrum of C2H4 at a column density of 1.010^^ cm ^, excitation temperatures of 200 (top), 
500 (middle) and 1000 K (bottom), and spectral resolving powers of 600 (left), 3000 (middle) and 50,000 (right). 



Jeanette E. Bast et al.: Exploring organic chemistry in planet-forming zones, Online Material p 5 



T3 

N 

a 

o 



0.99 



1.00 



: 

T = 200 K 


1 1 

R - 600 : 


1 

: T = 500 K 




: T = 1000 K 



0.99 



0.99 



0.97 



R ~ 3000 




1.0 



10 



12 



8 10 12 14 

Wavelength [/xm] 



0.8 
1.0 



0.9 



0.8 
1.0 



0.9 



0.8 



If 



R ~ 50,000 



nr 




10 



IS 



Fig. B.7. The synthetic spectrum of NH3 at a column density of 1.010^^ cm ^, excitation temperatures of 200 (top), 
500 (middle) and 1000 K (bottom), and spectral resolving powers of 600 (left), 3000 (middle) and 50,000 (right). 
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Fig. B.8. The synthetic spectrum of SO2 at a column density of 1.0' 10^^ cm~^, excitation temperatures of 200 (top), 
500 (middle) and 1000 K (bottom), and spectral resolving powers of 600 (left), 3000 (middle) and 50,000 (right). 
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Fig. B.9. The synthetic spectrum of C2H6 at a column density of 1.010^^ cm ^, excitation temperatures of 200 (top), 
500 (middle) and 1000 K (bottom), and spectral resolving powers of 600 (left), 3000 (middle) and 50,000 (right). 
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Fig. B.ll. The synthetic spectrum of CqHq at a column density of 1.010^^ cm ^, excitation temperatures of 200 (top), 
500 (middle) and 1000 K (bottom), and spectral resolving powers of 600 (left), 3000 (middle) and 50,000 (right). 
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Fig. B.12. The synthetic spectrum of C3H4 at a column density of l.OTO^^ cm ^, excitation temperatures of 200 (top), 
500 (middle) and 1000 K (bottom), and spectral resolving powers of 600 (left), 3000 (middle) and 50,000 (right). 
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Fig. B.13. The synthetic spectrum of HC3N at a column density of 1.010^^ cm ^, excitation temperatures of 200 (top), 
500 (middle) and 1000 K (bottom), and spectral resolving powers of 600 (left), 3000 (middle) and 50,000 (right). 
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Fig. B. 14. The synthetic spectrum of CH3CN at a column density of 1.0 10^^ cm~^, excitation temperatures of 200 
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